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ABSTRACT 


The  forces  and  moments  acting  on  a 
submersible  are  computed  when  it  is 
beneath  the  free  surface  or  near  a  wall. 
The  method  used  in  this  report  is  based 
on  potential  theory.  In  the  computation 
of  forces  and  moments  on  the  hull,  a 
three-dimensional  method  is  applied.  The 
free  surface  condition  is  linearized  and 
the  body  boundary  condition  is  exact. 

The  body  surface  is  discretized  with 
surface  elements  and  the  singularity  of 
source  and  sink  is  distributed  on  them. 
The  strengths  of  the  unkown  sources  and 
sinks  are  determined  through  the  body 
boundary  condition.  A  two-dimensional 
method  is  applied  compute  forces  and 
moments  of  the/control  surfaces.  The 
boundary  conditions  are'.same  as  those  for 
three-dimensional  case.  Computed  forces 
and  moments  of  the  control  planes  are 
added  to  those  of  the  bare  hull.  To 
include  the  interference  effect  of  the 
hull  on  the  control  planes,  the  flow 
velocity  at  each  control  plane  is 
computed  with  existence  of  the  hull. 

There  are  some  discrepancies  between 
computed  results  and  experimental  data 
because  of  the  effect  of  viscosity.  The 
overall  trends  in  the  computed  results 
are  same  as  those  of  the  experiments. 
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INTRODUCTION 


There  Is  a  need  to  develop  an  analytical  method  to  compute  the  forces  and 

moments  acting  on  a  submersible  when  it  is  beneath  the  free  surface  or  near  a 

wall.  The  methodologies  presently  available  to  compute  these  hydrodynamic 

forces  and  moments  are  limited  in  application.  Pond  *  computed  the  moments 

developed  on  a  Rankine  ovoid  by  using  the  method  of  axial  distribution  of 

2 

sources  and  sinks.  McCreight  improved  the  method  used  by  Pond  by  distributing 
the  dipoles  and  computed  the  vertical  force.  Their  results  are  in  good  agreement 
with  the  experimental  data.  However,  the  methods  used  by  Pond  and  McCreight  are 
valid  only  for  zero  pitch  angle. 

The  present  method  computes  the  vertical  and  longitudinal  forces  and  pitch 
moment  acting  on  a  submersible  with  and  without  control  planes  when  the 
submersible  moves  beneath  the  free  surface  or  near  a  wall.  This  method  is  more 
accurate  mathematically  than  that  of  Pond1.  Furthermore,  the  effect  of  pitch 
angle  is  incorporated  in  the  computation.  For  the  computation  of  forces  and 
moment  acting  on  the  bare  hull,  the  so-called  panel  method  is  applied.  The  body 
surface  is  discretized  with  many  quadrilateral  planes,  and  sources  and  sinks  are 
distributed  on  these  surface  elements.  This  method  can  be  applied  to  a 
submerged  body  of  arbitrary  shape.  To  compute  the  forces  and  moment  of  the 
control  planes,  the  method  developed  by  Giesing  and  Smith  is  used. 
Two-dimensional  sources  and  sinks  are  distributed  around  the  sections  of  a 
control  plane,  and  a  vortex  located  in  the  middle  of  each  section  is  introduced 
to  compute  lift  force.  The  inlet  velocity  to  the  control  planes  is  computed  at 
the  tips  of  the  sections  of  the  control  plane  to  include  the  interaction  effort 
between  the  bare  hull  and  control  planes. 


♦References  are  given  on  page  Ub. 


VELOCITY  POTENTIAL 


The  coordinate  system,  oxyz,  moves  at  a  speed  U,  which  is  the  mean  forward 
speed  of  the  submersible  ulong  the  positive  ox-axis  (see  Figure  1).  The  posi¬ 
tive  oz-axis  is  always  directed  vertically  upwards.  The  origin  0  is  located 
above  the  center  of  gravity  of  the  submersible.  The  oxy-plane  is  the  plane  of 
the  undisturbed  free  surface. 

We  assume  initially  that  the  submersible  and  coordinate  system  are  sta¬ 
tionary  and  that  the  fluid  around  the  submersible  moves  toward  the  negative  ox- 
axis  with  uniform  speed  U.  Then,  the  total  velocity  potential  for  the  bare  hull 
can  be  expressed  by 

<J>(x,y,z)  -  -Ux  +  <t>(x,y,z)  O) 

where  4>  is  disturbance  velocity  potential  due  to  the  submersible.  The  distur¬ 
bance  velocity  potential  satisfies  the  following  conditions: 

1.  Laplace  equation  in  the  fluid  domain 

2  2  2 

+  iLi  +  L±  .  o  (2) 

2  2  2 
dx  dy  dz 

2.  The  linearized  free-surfacc  condition 

^4+ko^-0  (3) 


3. 


The  body  boundary  condition 


4.  The  radiation  condition:  the  disturbance  vanishes  sufficiently  fast 
far  ahead  of  the  ship. 


—  2  2  2 

im  R(  $  +  4  +  $  ) 

x  y  vz 

R  -  (x2  +  y2)1/2 


5.  The  bottom  condition:  the  normal  velocity  at  the  bottom  is  zero. 


,(x,y,z)-0  asz 


In  equation  (3),  the  constant  is  given  as 


k  -  & 

kO  2 

U 


and  in  Equation  (4),  the  unit  vector  which  is  oriented  normal  to  the  fluid  is 
given  by 


(n^ ,  ^2*  n3^  **  ^ 


<n4»  n5.  n6>  =*  rG  *  « 


where  r^  is  the  position  vector  of  the  center  of  gravity. 


The  solution  of  Equation  (2)  is  given  by  Brard  in  an  integral  form  as 


4>(x,y,z)  -  -  ~  If  G( P  ,Q)o(Q)dS(Q) 
4n  s 


where  P(x,y,z)  is  the  field  point,  Q(x0,y0,z0)  source  point,  S  the  wetted 

surface  of  the  body,  a  the  unknown  strength  of  sources  and  sinks  distributed  on 

the  body  surface,  and  G  is  the  Green  function  which  is  given  by  Wehausen  and 

,  ,  5 

Laltone  as 


G( P ,Q )  -  —  -  — - -I  sec  9  d9  / 

r  r  j  it  -n  o 

where 

r2  -  (x  -  xQ)2  +  (y  -  y0)2  +  (z  - 

2  2  2 

r  ”  (x  -  *0>  +  (y  -  y0)  +  (z  + 


-(z+z0)u  iliJU 

- - - - du 

2 

(10 

u  -  k0sec  9 

zo>2 

(11 

zo>2 

(12 

and 

(jj  ■  (x  -  x0)  c os 9  +  (y  -  yQ)  sin©  (13: 

Only  the  real  part  of  Equation  (10)  Is  used  in  later  computation. 

The  strength  of  the  sources  and  sinks,  a,  in  Equation  (9)  can  be  found  by 
substitution  of  Equation  (9)  into  Equation  (4)  to  obtain 


4x 


G(P,Q)o(Q)ds(Q) 


-Un, 


(14; 


The  solution  of  Equation  (14)  is  only  feasible  with  the  help  of  a  numerical 
procedure  which  will  be  given  later. 

The  force  acting  on  the  body  is  expressed  as 


"// 

S 


p  n  ds 


(15) 


where  p  is  the  pressure  around  the  body  and  is  linearized  from  Bernoulli's 
equation  as 


5 


(16) 


P  -  P  U  ♦* 

where  p  is  the  density  of  the  fluid. 

CONTROL  PLANE 

The  velocity  potential  for  the  control  plane,  whose  coordinate  system  is 
shown  in  Figure  2,  is  expressed  in  the  two-dimensional  domain  as 

4<(x,z)  »  -Ux  +  4>(x,z)  (17; 

where  <J >  is  the  disturbance  velocity  potential  due  to  a  section  of  the  control 
plane.  This  disturbance  potential  satisfies  the  following  conditions. 


Laplace  equation  in  the  fluid  domain 

+  lL±  .  o 
2  2 
dx  dz 


2.  The  linearized  free-surface  condition 


+  k  £4.o 

o  kO  u 

dx  dz 


3.  The  body  boundary  condition 


£4  - 


4.  Kutta  condition:  the  velocities  at  the  trailing  edge  elements,  one  at 
the  top  of  the  surface  and  one  at  the  bottom,  are  equal. 


In  addition,  the  disturbance  potential  should  satisfy  radiation  and  bottom 
conditions  similar  to  Equations  (5)  and  (6). 
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The  solution  of  Equation  (18)  Is  given  by  Gleslng  and  Smith  •  The  velocity 
potential  function  Is  divided  into  three  parts  as 


4>(x,z)  -  <Jij  +  r(4>2  +  4»3). 


The  velocity  potential  Is  due  to  sources  and  sinks  which  are  distributed  on 
each  section  of  the  control  plane.  <\>2  Is  the  velocity  potential  due  to  a  vortex 
of  unit  strength  which  Is  located  at  the  inside  of  a  section.  The  vortex  Is 
Introduced  to  calculate  the  lift  force  around  the  section.  4^,  the  last  term  of 
Equation  (21)  Is  the  velocity  potential  due  to  sources  and  sinks  of  unit 
strength  distributed  on  the  section.  This  velocity  potential  cancels  the  normal 
velocity  generated  by  4>2«  If  we  substitute  Equation  (21)  into  Equation  (20), 
the  following  conditions  are  given 


These  three  velocity  potentials  are  expressed  as 


<l>t(P)  *  /  °i(q)G2^P  *q)dl(q)  for  1*1  and  3 


<K(P)  =  Re{  i-  ln[x-a  +  i(z-b)]  -  4—  In  [(x-a)  +  l(z+b)] 


-2pv  / 


-ik[ x-a  +  i( z+b)  ] 


k  -  k. 


+  2x1  e'lkofx"a  +  lU+t>)]j 
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where  cf±  is  the  two-dimensional  strength  of  sources  and  sinks,  (a,b)  the 
location  of  the  vortex,  q(x0,zQ)  a  two-dimensional  source  point,  and  G 2  is  the 
two-dimensional  Green  function  given  by  Wehausen  and  Laitone5  as 


G2(p»q) 


Re  {in  [x-x0  +  i(z-z0)]  +  ln[x-xQ  +  i(z+zQ)] 

(26) 

«  e-i^[x-x0  +  i(z+zQ)]  dk  _  2itl  e-ik0[x-x0  +  i(z+zQ)]l 

+  2pv  / 

o  ° 


In  Equations  (25)  and  (26),  Re  denotes  the  real  part  of  a  complex  quantity  and 

pv  the  principal-value  integral.  The  vortex  strength  T  in  Equation  (21)  will  be 

determined  with  the  Kutta  condition.  Tne  pressure  around  each  section  of  the 

control  plane  can  be  obtained  from  Equation  (16)  by  substituting  ib  for  a  The 

x  x 

force  and  moment  acting  on  the  section  can  be  computed  with  an  equation  of  the 
same  form  as  Equation  (15). 
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Si-  o 

bz 


at  z  ■  0 


(29) 


All  other  boundary  conditions  are  the  same  as  those  described  in  Equations  (2), 
(4)  -  (6),  (18)  and  (20).  The  solutions  in  this  case  have  the  same  form  as 
Equations  (9)  and  (21),  with  different  Green  functions.  The  Green  functions  can 
be  expressed  by  application  of  the  method  of  images  (see  Reference  6)  as 


«r.Q>  "  7  *  7, 


(30) 


for  the  case  of  bare  hull  and 


G,(p,q)  "  Re{ln[x-x0  +  i(z-zQ)]  +  In  [x-x0  +  i(z+zQ)]} 


(31) 


for  the  case  of  control  plane.  The  velocity  potential  due  to  a  unit  vortex 
located  at  (a,b)  should  be 


<j,_(p)  *  Re{—  ln[(x-a)  +  i(z-b)]  -  In  [(x-a)  +  i(z+b)]} 


2* 


2x 


(32) 


To  find  the  velocity  potentials  for  the  wall  condition,  Equations  (10),  (25)  and 
(26)  should  be  replaced  by  Equations  (30),  (32)  and  (31),  respectively. 
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NUMERICAL  PROCEDURE 

To  determine  the  unknown  strength  of  the  sources  and  sinks  in  Equation  (14), 
we  first  discretize  the  wetted  surface  of  the  body  S  with  many  quadrilateral 
elements.  Figure  3  shows  the  discretized  surfaces  of  a  spheroid.  Furthermore, 
we  assume  that  the  unknown  strength  of  the  sources  and  sinks,  o(Q),  is  constant 
on  each  surface  element.  Then  Equation  (14)  can  be  written  as  follows 


N 


h  jij  h  G(pi>QJ)  ssj  '  -""u 


for  i  *  1,  2, 


(33) 


Once  the  derivations  of  the  Green  function,  G,  are  evaluated  numerically, 
Equation  (33)  can  easily  be  solved  for  Oj.  The  numerical  evaluation  of  the 
Green  function  and  its  derivatives  is  the  most  difficult  part,  and  consumes  a 
lot  of  computer  time  when  running  the  program. 

The  normal  derivative  of  the  Green  function  in  Equation  (33)  can  be 
expressed  as 


bG 

bn 


bG 

bx 


bG  bG 

n,  +  —  n0  +  —  n, 
1  by  ^  bz  J 


(34) 


The  task  is  to  find  three  derivatives  of  the  Green  function.  To  do  this,  we 
first  derive  the  Green  function  suitable  to  numerical  evaluation.  The  Green 
function.  Equation  (10),  can  be  rewritten  as 


G  “7 


—  _£  [ 
r!~  11  o 


sec  9  d8  {/ 
o 


-zu,  iui.u  ,  i£5_u) 

e  (e  +  +  e 


u  -  u. 


du 


(Lj) 


CD  -zu,  -iw+u  ,  -l£>_u , 

+  J  *  ■+-  « - 1  du} 

O  u  ■  U° 

(L2) 


(35) 
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where 


“+  ■  <*-xo)  co8  0  ♦  <-y-yQ) 


sin  6  ■  u> 


5_  »  (x-x0)  cos  8  -  (y-yQ)  sin  8 


(36) 


z  *  z  +  z. 


Jo 


The  notation  under  two  of  the  Integrals  In  Equation  (35),  namely,  or  L^, 

represents  the  integral  paths  shown  in  Figure  4.  To  derive  the  derivatives  of 

the  first  two  terms  of  Equation  (35),  the  method  which  Hess  and  Smith  developed 

8 

for  unbounded  flow  can  be  used,  Hong  and  Paulling  have  applied  this  method  for 
the  computation  of  motions  of  a  body  in  waves,  and  their  paper  can  be  examined 
for  the  details.  Since  Reference  8  can  be  examined  for  these  details,  we  will 
present  here  the  derivation  of  the  last  two  terms  of  Equation  (35). 


We  write  Equation  (35)  as  follows 


G  -  -  -  -  +  G, 
r  rt  3 


(37) 


Chen,  Che  InCegral  pach  becomes  that  shown  in  Figure  5.  The  integral  path, 
la  now  transferred  onto  Che  pach  C,  in  Figure  5  as 


/  ' 
o 

(Lx) 


-ZU  l<i)_u  ~v  -V 

e  6  du  -  e"uo(*  -  lw>  J  £— 

'  tt 


U  -  Uf 


dv 


Cl 


(39) 


The  integral  over  Cj  can  be  easily  evaluated  with  the  method  of  residue  as 


(40) 


with  the  introduction  of  the  exponential  Integral 


the  integral  over  Cj  becomes 


-v 


u  <  o 
U)  >  o 


(42) 
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By  substituting  Equation  (37)  Into  Equation  (33),  the  Integral  over  Lj  is  given 

by 


—  <v 

-ZU/  ilO+U  .  1w_Uv  V 

re  (e  +  +  e  )  .,  -urt(z  ~  lu>+)  t„  r 

i  - -  du  -  e  o  +  {Ej  [-u0(z  -  lw+)J  - 


+  e-uo(*  ~  lw-> 


(Ejt-u^r  -  is.)]  ♦  {2,1|,  "+•  <  ° 

0  u»+,  ti>«  >  o 


For  the  Integral  path,  L_,  we  change  the  variable  as  follows: 


v  -  (u  -  uQ)  (z  +  la)  with  u>  -  w+  or  u>_; 


then,  the  new  Integral  path  Is  given  In  Figure  6.  With  this  new  variable,  the 
integral  path,  Lj,  is  now  C^.  Employing  the  same  procedure  as  used  previously 
to  evaluate  the  integral  over  Lj,  the  integral  over  becomes 


•v  »v 

•  -Tu,  lu+u  .  -lu_ux 

I  £  is-  -  ■»  e - 1  dU 


«“o<*  *  iS*>  is+»i  *  2;‘i 


♦  e-uo'z  *  ‘"-Vjt-UoU  ♦  lS_)]  ♦  2,11,  <  ° 

0  3+,  w_  >  o 


w_  <  o  (44) 


Iv.vV-v 


Substitution  of  Equation  (43)  and  (44)  into  Equation  (38),  allows  us  finally  to 
express  in  a  convenient  form  for  numerical  evaluation  as 


n/2 


G3“ 


“4  I  uQ  d0  {e  U°^Z  iW+){E1[-u0(z  -  iw+)]  -  2lti} 

0 


J 

0 


+  e_uo(z  ”  la>-) 


{E  [uQ(z  -  iuL)]  -  2nl} 

0 


+  e-u0(z  -  iu+)  {Ei[^o(7+  lS+)]  +  2^i} 


+  e-uQ(z  +  iw_)  (Ei[^o<7+  lSL)]  +  2ni}}>  y  Z.  <  0 

0  u>+,  u>_  >  0 


If  we  let 


uQe  U°^Z  iW+^  Ej[-uq(z  -  iu>+)]  -  R  +  il 


u0e"Uo(z  '  lu,+  )  -  C  ♦  iS 


u„e 


-UQ(Z  -  i(i)_) 


E^[u0(z  -  iw_) ]  =  R  +  il 


u0e 


-u0(z  -  iu>_)  _  c  + 


C  +  IS 
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(45) 


(46) 


Si 


i 


/ 

A 


'i 


■Ji 


we  can  express  as 


g3  -  -  f  f  [R  +  (2J)S  +  R  +  (2J)s]d0f  tJ+’  W‘  <  ° 


G+,  G-  >  0 


The  derivatives  of  G^  can  be  easily  expressed  from  Equation  (38) 


where 


1  u  cos 9  G, 


■  1  u  sinG  G.  8ign(k)) 


3G3 

-  ■  u  G- 

bz  3 


sign  (u)  *  1  for  u>  -  u>+ 


0)  *  u_ 


With  substitution  of  Equation  (45)  into  Equations  (48),  (49)  and  (50),  the  deri¬ 
vatives  of  Gj  can  be  expressed  in  final  form  as 


-2  ^  ~2 
\  +  “+ 


i  +  ( 2u)c 

0 


+  (J~  -  I  +  (27t)  C}d9 

-2  ~2  0 
z+  +  o)_ 


-  I  j  8lne  t-  ~ 


-2  ,  ~2 
z  +  w 
+  + 


I  +  (2n)C 
0 


+  M-  +  I  -  (2lX)  c}d0 
-2  ~2  0 
z  .  +  0J 

*T  “ 


3  2 

ST"'*  J  u° 


- 5 -  +  R  +  (2^)S 

_2  ~2  0 
z  +  w 
+  + 


- - +  r  +  (2*)  side 


^rwwinr flrwri  fn" 


WMM  R! I WW IW W* U *  7* UTT*  V*\Pf  V* Vt  V¥V*  \rv  \!*  ^  l*  wvtvjtwvt  l«sr  irmr «n  v 


In  these  equations,  the  upper  value  is  taken  when  Q+  or  Q_  is  larger  than  and 

equal  to  0  and  the  lower  value  is  taken  when  Q.  or  0  is  negative.  Inglis  and 

9  + 

Price  have  derived  the  Green  function  for  a  translating  and  pulsating  source 

s'ftilar  to  Gj  with  the  change  of  variable  and  with  the  change  of  integral  path. 

When  the  argument,  -uQ(z+iu)  or  -u0(z-iQ>)  is  small,  the  exponential  integral 
can  be  expressed  in  a  series  form  (see  Reference  10)  as 

E.(z)  -  -y  -  In  z  -  l  z°.  |  arg  z|  <  x  (55) 

1  n-1  nn! 


where  z  is  any  complex  number  and  y  m  0.5772157  is  Euler's  constant.  When  the 
argument  becomes  very  large,  the  exponential  integral  is  evaluated  by  the  method 
of  Todd.  For  a  large  complex  number  z  -  x  +  iy,  the  exponential  integral 
multiplied  with  exponential  function  is  expressed  as 


e  Ej ( z )  -  II  -  il2 


where 


a> 

I  -U  X  +  U 

I  6  2  2 

J  (x  +  u)  +  y 


du 


X2  ” 


QD 

/•"TTr^r 


du 


(56) 


(57) 


(58) 
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For  the  evaluation  of  Ij  and  I2»  we  use  the  following  approximation 


I 


J  e_t  f(t)  dt  H  I  X<n)  f(x^n)) 
0 


(59) 


where  x^°^  are  the  zeros  of  the  Laguerre  polynomial  and  X^n^  are  the 

corresponding  Christoffel  numbers.  In  this  report,  n  is  taken  as  5,  and  x. 

(5)  1 

and  X  are  given  in  Table  I. 


(5) 


Table  I: 

Zeros  of  Laguerre  Polynomials  and 
Christoffel  Numbers 


0.2635603 

1.4134031 

3.5964258 

7.0858100 

12.6408008 


0.5217556 

1.3986668 

0.0759424 

0.0036117 

0.0000234 


The  exponential  functions  in  Equations  (47),  (52),  (53)  and  (54)  are  highly 
oscillatory  when  the  real  part  of  the  complex  argument  Is  small  and  the  imagi¬ 
nary  part  is  large.  One  of  complex  arguments  in  Equation  (46)  can  be  expressed 


with 


R2  -  (x  -  xQ)2  +  (y  -  y0)2 


0  -1  ,1  ~  Yov 

P  ■  tan  (— - ■) 

x  -  xQ 


When  ©2  -  9j  is  small,  the  Integral  between  Qj  and  ©2  can  be  approximated  as 


0 


J 


2 

uQ(C  +  iS)d0 


’1 


w 

; 

9i 

e 

•f 


kQ  sec  9 


exp{-kQ  sec20[z  +  1R  cos(0  +  0)]}  (61) 


(A  +  IB)  exp[ l(a0  +  b)]d0 


where  A,  B,  a,  and  b  are  constant  between  0^  and  02*  The  last  expression  In 
Equation  (61)  can  be  analytically  evaluated.  Once  the  normal  derivative  of  the 
Green  function,  which  contains  Equation  (61),  is  numerically  evaluated,  Equation 
(33)  can  be  expressed  as  a  system  of  linear  equations  and  can  be  solved  by 
the  method  of  Gauss  elimination. 


For  the  computation  of  forces  and  moment  of  a  control  plane,  we  discretize 
each  section  of  the  control  plane  with  straight  segments  as  shown  in  Figure  7. 
The  unknown  strength  of  the  two-dimensional  sources  and  sinks  on  each  segment  is 
assumed  to  be  constant.  The  subscript  l  for  in  Equation  (27)  will  be  dropped 
from  now  on  to  avoid  confusion  associated  with  that  of  segments.  a  represents 
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each  of  and  In  the  following  equations*  Equation  (27)  can  be  approximated 


i  M  C  a 

2*  j-i  °(qJ)  J  ^  G2  (Pi’  V  dSJ(qJ>  '  U  nii’ 


for  1  ■  l ,  2 , . . .  ,M 


the  normal  derivative  of  the  two-dimensional  Green  function  is 


2>G2  SG2  $G2 

ftn  dx  nl  +  dz  n3 


We  use  the  complex  variables  in  derivation  of  derivatives  of  the  Green  function 
suitable  for  numerical  evaluation.  If  we  introduce  the  following  complex 
variables 

C  *  x  +  12  and  C  -  xQ  +  lz0,  (64) 


G2(p,q)  -  Re{F(C,C0)} 


where 


F(C,C0)  =*  ln( C  -  CQ)  +  ln( C  -  C0)  +  pv 


-  2x1  e“ik°(C  "  Cq) 


In  Equation  (66),  C  is  the  complex  conjugate  of  C  . 

o  o 

in  Equation  (63)  can  be  expressed  as  functions  of  F 
ac2  n  ,&F, 


e-ik(C  -  CQ) 
k  -  kn 


Then,  the  two  derivatives 


JL.V,  Ox.',  OO  r~.  O  «  ^  a  *  -  -  O  V .  -  «  ' *  *  A*  »  I*/1.  -  w  *  a.y« 


The  derivative  of  Equation  (66)  with  respect  to  C  is 


GD 

9F  _  1  ,1  i  f 

—  -  - —  + - pv 

ac  C  -  Co  C  -  Co  x  I 


-ik(C  -  CQ) 


k  -  k. 


-  2n  kQ  e 


-i  kn(C  -  Crt) 


Integrating  the  last  equation  with  respect  to  sj (qj )  we  obtain 


/  ~  dsj(q-j)  -  -(cos  aj  -  i  sin  aj)[ln  r^+1  +  i0j+1  -  In  r^  -  19^] 


-  (cos  aj  +  1  sln  aj)tln  rj+i  +  i0j+l  ~  ln  rj  "  l6j 1 


{  (cos  Qj  +  i  sin  V(PvJ^ 


k(z  +  z_  ) 

dk  °j+l  r  , ,  i 

e  Lc°8  k(x  “  xoj+i) 


,  f  dk  k(z  +  zoi\ 

i  sin  k(x  -  x0^+1)J  -  pv  |  — - —  e  J  [cos  k(x  -  x0^ ) 


* 


niiv.uyp  mwyrmymJmyyy  ’r^K’/.v.v.v.v; v.  ^ 


^1^ 


1TWYWT 


-  i  sin  k(x  -  xQ j ) ] } 

+  2xi  (cos  <xj  +  i  sin  aj){e  °J+1^[cos  kQ(x  -  xQj+j) 

-  1  sin  kQ(x  -  x0  )] 

-  e  °  °J+1  [cos  kQ(x  -  x0^)  -  i  gin  kQ(x  -  xQj )  ] }  (70) 

the  new  variable  and  notation  are  given  in  Figure  8.  The  details  for  derivation 
of  Equation  (70)  are  given  by  References  3  and  12. 

The  numerical  procedure  for  determining  and  T  is  as  follows:  First,  we 
evaluate  the  normal  derivative  of  the  two-dimensional  Green  function  which 
enables  us  to  solve  Equation  (62)  for  o^.  Next,  we  calculate  the  normal 
velocities  at  segments  due  to  a  unit  vortex  located  at  the  center  of  the 
section.  To  eliminate  these  normal  velocities,  we  distribute  a ^  at  each  segment 
of  the  section  and  determine  w*th  Equation  (23).  We  finally  solve  for  the 
vortex  strength  T  in  Equation  (21)  using  the  Kutta  condition. 

The  computational  procedure  of  forces  and  moments  is  as  follows:  The  forces 
and  moments  of  the  bare  hull  are  first  computed.  This  means  it  is  assumed  no 
interference  of  control  planes  to  the  bare  hull.  To  Include  the  interference 
effect  of  the  bare  hull  to  the  control  planes,  the  flow  velocity,  U,  in  Figure  2 
is  different  from  U  in  Figure  1.  A  control  plane  is  cut  at  three  different 
spanwise  locations.  At  the  leading  edge  of  each  location  (or  section),  U  is 
computed.  This  velocity  is  different  from  the  steady  forward  speed,  U,  in 
Figure  1.  Two-dimensional  forces  and  moments  are  computed  at  each  section;  and 
these  are  summed  along  the  spanwise  direction.  Finally,  these  summed  forces  and 
moments  are  added  to  those  of  the  bare  hull. 
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NUMERICAL  RESULTS 


In  this  section  we  will  calculate  the  forces  and  moments  acting  on  a  submer¬ 
sible  moving  under  the  free  surface  or  under  ice.  Numerical  results  based  on 
the  present  method  will  be  compared  with  results  obtained  from  earlier  theoreti¬ 
cal  methods  and  experiments.  First,  forces  acting  Independently  on  the  bare 
hull  and  control  plane  will  be  compared;  then,  results  for  the  combined  hull  and 
plane  will  be  compared. 

Figure  9  shows  the  forces  and  moments  on  a  Ranklne  ovoid  having  a  length  of 
4  ft  (1.22ra).  The  ratio  of  length  to  diameter  (D)  is  10.5.  The  Ranklne  ovoid 
is  created  by  distributing  sources  and  sinks  along  a  line  which  Is  assumed  to  be 
In  a  uniform  stream  In  an  unbounded  fluid  (see  Reference  12).  The  computed  ver¬ 
tical  forces  ( F ' )  are  generally  somewhat  larger  than  the  experimental  data  for 
z 

smaller  submergence  (h)  and  smaller  for  larger  submergence.  The  computed 
moments  (M' )  are  generally  smaller  than  the  experimental  data  for  all  sub- 

y 

mergences  when  the  Froude  number  la  larger  than  0.5.  However,  the  overall  trend 
of  the  computed  results  Is  similar  to  that  of  the  measurements.  Results  calcu¬ 
lated  by  the  present  method  are  in  good  agreement  with  values  obtained  from  pre¬ 
vious  analytical  methods.  Figure  10  shows  the  same  results  plotted  as  a 
function  of  submergence.  As  the  submergence  becomes  larger  than  three  times  the 
diameter,  the  force  and  moments  decrease  rapidly  for  all  Froude  numbers. 

Figures  11,  12  and  13  show  similar  data  for  a  spheroid  whose  ratio  of  length 
to  diameter  is  7.  The  findings  here  are  somewhat  different  from  those  discussed 
previously  for  the  Ranklne  ovoid  (see  Figure  9).  When  pitch  angle  is  zero  (see 
Figure  11)  the  experimental  data  for  the  spheroid  extracted  from  Reference  13 
are  larger  than  the  computed  values  at  Froude  numbers  less  than  0.5.  The 
agreement  between  computation  and  experiment  is  generally  good  when  the 
submergence  is  equal  to  or  larger  than  the  diameter:  when  the  ratio  of 
submergence  to  diameter  Is  0.75,  there  is  substantial  difference.  When  the 


plcch  angle  ie  2.5  degrees  (bow  up),  the  experimental  results  are  larger  than 
those  from  computation  for  essentially  all  Froude  numbers  (see  Figure  12).  As 
shown  in  Figure  13,  negative  pitch  angle  of  2.5  degrees  (bow  down),  the 
prediction  and  measurement  of  vertical  force  are,  for  the  most  part,  in  better 
agreement  than  those  for  zero  or  positive  pitch  angle.  The  computed  moments  are 
larger  than  the  experimental  results  for  all  submergences  and  Froude  numbers. 

These  discrepancies  for  non-zero  pitch  angle  might  be  caused  by  the  viscosity  of 
the  fluid  and  also  vortex  shedding  from  the  body  or  the  cross  flow  effect. 

Figure  14  shows  the  results  of  lift  for  a  hydrofoil  whose  section  has  the 
shape  of  NACA  4412.  The  results  shown  in  this  figure  are  those  for  a  two- 
dimensional  section.  The  agreement  between  computation  based  on  the  present 
method  and  experimental  data  from  Reference  14  is  good  when  the  angle  of  attack 
i 8  smaller  than  4  degrees  (see  bottom  figure).  As  the  angle  of  attack 
increases,  the  computed  results  become  larger  than  the  measurements  particularly 
at  the  lower  Froude  number.  It  is  of  interest  to  observe  that  for  this  sub¬ 
mergence  the  lift  coefficient  at  the  higher  Froude  number  is  smaller  than  that 
at  the  lower  Froude  number.  When  the  angle  of  attack  is  10  degrees  and  the 
Froude  number  Is  0.922  (see  top  figure)  there  is  substantial  discrepancy  between 
the  computed  and  experimental  results  as  the  submergence  increases.  The  results 
computed  by  the  present  method  are  further  compared  with  the  method  developed  by 
Wadlin  and  Christopher  In  Reference  17,  Equation  (15).  The  calculations  were 
made  for  an  aspect  ratio  of  10  to  stimulate  a  two-dimensional  lift  coefficient. 

It  should  be  pointed  out  that  Equation  (15)  of  Reference  17  is  strictly 
applicable  for  aspect  ratio  of  0.125  to  10.  Their  results  are  almost  half  of 
those  by  the  present  method.  Compared  with  the  experiment,  the  results  of 
Wadlin  and  Christopher  show  better  agreement  than  those  by  the  present  method. 

Figure  15  shows  the  results  of  lift  for  a  hydrofoil  whose  section  has  the 
shape  of  NACA  64A010.  The  aspect  ratio  of  this  hydrofoil  is  6.  As  shown  in  the 
bottom  figure,  for  small  angle  of  attack  there  are  small  differences  between 
computation  and  the  experimental  results  taken  from  Reference  15.  The  differen¬ 
ces  increase  as  the  angle  of  attack  increases.  The  reason  for  these  discrepancies 
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may  be  that  first,  the  effect  of  viscosity  is  not  included  in  the  computation, 
and  second,  the  test  was  done  with  turbulence  stimulation.  The  top  figure 
indicates  that  for  an  angle  of  attack  of  4  degrees,  the  computed  lift 
coefficient  is  about  twice  the  measured  value  when  the  submergence  is  large. 

When  the  ratio  of  submergence  to  chord  is  smaller  than  0.5,  the  agreement 
between  computation  and  experiment  is  good.  For  small  submergences,  the  lift  is 
more  affected  by  the  free  surface  than  viscosity,  and  when  the  submergence  is 
large,  the  reverse  is  true.  On  the  other  hand,  compared  with  the  experiment, 
the  results  of  Wadlin  and  Christopher  show  better  agreement  than  those  by  the 
present  method  when  the  submergence  ratio  is  larger  than  2. 

When  the  location  of  the  sail  plane  or  stern  plane  of  a  submersible  is  at  a 
depth  of  two  or  three  times  of  chord,  a  better  lift  computation  can  be  expected 
with  the  method  developed  in  Reference  17. 

Figure  16  shows  forces  and  moment  for  Model  4621  with  and  without  stern- 
planes  at  deep  submergence.  The  computed  and  experimental  axial  forces  are 
fairly  steady  at  different  trim  angles.  The  vertical  forces  for  the  bare  hull 
are  computed  to  be  significantly  less  than  those  of  experiment.  However,  the 
results  of  vertical  forces  with  sternplane  agree  very  well  with  those  of  experi¬ 
ment  except  at  <x=12°.  The  reason  for  the  good  agreement  for  the  case  with 
sternplane  is  that  the  vertical  force  (in  this  case  lift)  of  the  control  plane 
alone  is  over-estimated  at  deep  submergence  as  shown  in  Figures  14  and  15;  and 
this  over-estimation  is  compensated  with  the  under-estimation  of  the  vertical 
forces  of  base  hull.  The  moments  are  computed  to  be  larger  than  the  experimen¬ 
tal  values. 

Figure  17  shows  the  results  of  vertical  force  of  a  spheroid  with  L/D=7  near 

1  8 

a  wall.  The  computed  results  are  compared  with  the  results  of  Newman. 

Newman's  method  was  developed  using  slender  body  theory  and  with  the  assumption 
of  L/D  >>  1  and  h/D  <<  1.  The  vertical  forces  computed  by  the  present  method 
are  smaller  than  those  computed  by  Newman.  It  is  unkown  which  method  is  more 
accurate. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  forces  and  moment  acting  on  a  submersible  moving  beneath  the  free  sur¬ 
face  or  near  a  wall  with  a  pitch  angle  for  the  bare  hull  or  a  control  surface 
deflection  have  been  computed  using  potential  flow  theory.  There  is  good 
agreement  between  calculation  and  experiment  for  the  bare  hull  when  the  ratio  of 
submergence  depth  to  the  diameter  is  larger  than  1  and  the  pitch  angle  is  zero. 
For  the  control  planes  alone,  the  good  agreement  between  calculation  and 
experiment  is  achieved  when  the  ratio  of  submergence  depth  to  chord  is  small. 
When  this  ratio  is  large,  the  force  due  to  viscosity  is  dominant  and  the  present 
method  generally  computes  forces  larger  than  those  of  experiment. 

From  the  present  stuuy,  the  following  conclusions  may  be  drawn: 

1.  The  numerical  evaluation  of  the  Green  function  needs  to  be  improved 
when  the  ratio  of  submergence  depth  to  diameter  is  smaller  than  1.0. 
Near  the  free  surface,  the  Green  function  has  oscillatory  behavior  and 
the  exponential  function  decays  slowly. 

2.  For  nonzero  angle  of  attack,  there  is  a  cross  flow  or  vortex  shedding. 
The  effect  of  cross  flow  or  vortex  shedding  should  be  included  in  the 
computation  with  the  help  of  viscous  flow  analysis. 

3.  When  the  ratio  of  submergence  depth  to  the  mean  chord  of  a  control 
plane  is  large,  the  lift  is  dominated  by  the  contribution  of  viscosity. 
The  viscous  effect  on  the  lift  should  be  incorporated  in  the  future 
computat ion. 

A.  To  improve  the  results  of  potential  theory,  an  empirical  or  analytical 
approach  should  be  developed  with  inclusion  of  viscous  effect. 


ACKNOWLEDGEMENT 


The  author  acknowledges  the  support  of  Messrs.  M.  Davis  and  D.  Milne.  He 
also  expresses  his  thanks  to  Messrs.  J.  Feldman  and  T.  Moran  for  their  coor¬ 
dination  of  this  work.  Finally,  the  author  thanks  Mr.  A.  Gersten  for  his  exten¬ 
sive  editorial  advice. 


surface 


when  65  <0 


•u  (Z-iW) 
o 


■u  (Z-iG3) 

o 


when  U  >  0 


Figure  5  -  Change  of  Integral  Path  for  1 
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Figure  9  -  Force  and  Moment  on  a  Rankine  Ovoid  with  I  / 1)  =  10,') 
at  Different  Submergences 

30 


V.AV, 


Figure  14  -  Lift  of  a  Hydrofoil  of  NACA  441f’  Shape 
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Figure  15  -  Lift  of  a  Hydrofoil  of  NACA  64A010  Shape 
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